Accurate Predictions for the Scale-Dependent Galaxy Bias from Primordial 

Non-Gaussianity 
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The large-scale clustering of galaxies can serve as a probe of primordial non-Gaussianity in the 
Universe competitive with the anisotropics of the CMB. Here, we present improved theoretical pre- 
dictions which include an important, previously overlooked correction to the bias. We demonstrate 
that the new predictions are able to reproduce the results of N-body simulations, explaining the 
significant departures seen from previous theoretical results. These refined predictions open the way 
to accurate constraints on primordial physics with large-scale structure surveys. 
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Measurements of the temperature and polarization 
anisotropies of the cosmic microwave background (CMB) 
radiation strongly support inflation as the mechanism 
generating the initial curvature perturbations [1], here 
parametrized by the Bardecn potential perturbations 
(p = [2]. While the simplest, single-field slow-roll 
models predict nearly Gaussian fluctuations, many other 
inflationary scenarios yield a potentially detectable level 
of primordial non-Gaussianity, which can generally be de- 
scribed in terms of the higher-order {N > 3) correlation 
functions of 4> in Fourier space, C^^^- Because different 
seeding mechanisms of the initial curvature perturbations 
leave distinct signatures in the shape and amplitude of 
these higher-order correlation functions, a precise mea- 
surement of non-Gaussianity would strongly constrain 
the physics of inflation. 

As the primordial curvature perturbations eventually 
grow into the large scale structure observed today, their 
higher order correlation functions can be determined 
through a measurement of the clustering of dark matter 
halos [3] which host the galaxies. For example, primor- 
dial non-Gaussianity of the local /NL0^-type induces a 
strongly scale-dependent bias Abi(fc, z) in the clustering 
of halos [4] . The power spectrum Ph of halos of mass M 
at redshift z is related to the matter power spectrum Pm 
through 

Pnik, M, z) = [bf (fc, M, z)] z) (1) 

= [bf{M,z) + Ahi{k,M,z)Yp,n{k,z) 
2/nl'^c 



Abi(fc,Af,z) = (6f(M,z) - 1) 



M{k,z) 



(2) 



Here, bf{M,z) is the linear, Gaussian (Eulerian) bias of 
halos (6i = 6f — 1 is the corresponding Lagrangian bias), 
6c — 1.69 is the critical density for spherical collapse, and 
the function Mik,z) = 2k^T{k)D{z)/{3nmH§) depends 
on the linear growth factor D{z), the matter transfer 
function T{k), the present-day matter density parameter 
rim and the Hubble constant Hq. In what follows, we 



shall omit the mass and redshift dependence for concise- 
ness. 

The scale-dependent non-Gaussian (NG) bias expres- 
sion [Eq. (2)] may be extended (for = 3) using a 
peak-background split (PBS) argument [7]. Alterna- 
tively, Eq. (2) may be generalized to any model of primor- 
dial NG by computing the correlation function of high 
threshold regions [5]. Even though the result is strictly 
valid in the high-peak limit only, i.e. v = Sc/ctqs ^ 1 
where ctos is the variance of the density field smoothed 
on scale Rs{M) = (3M/47rp)^/^, it is usually extended 
to arbitrary peak heights upon identifying v/uqs with bi. 
For a given arbitrary N , this gives 



Ab{'^P'(fc) = 



(A^-1)! \ao 



N- 
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(k) 



Ms{k) 
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We have defined Ms{k) = M{k)WR^{k), where WR^{k) 
is the Fourier transform of the spherical tophat filter with 
radius Rs- Furthermore, we have introduced the shape 
factor 



J-W(fc) 
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with 5 = — fci — • • • — — k, and P,f,{k) being the power 
spectrum of the Bardecn potential 0. For the constant 
/nl0^ model, J-i^\k)/ f^i^ is equivalent to the function 
J-R{k) defined in [5]. As Fs (k) — >■ /nl on large scales, 
Eq. (3) reproduces Eq. (2) in the limit — >■ 0. Since, in 
the limit ^ 1, both thresholding and PBS approaches 
yield Eq. (3), we call it the high-peak result. 

Measurements of the clustering of dark matter ha- 
los in non-Gaussian N-body simulations have confirmed 
the vahdity of Eq. (1). The fc-dependent NG bias cur- 
rently yields constraints on the nonlinear parameter /nl 
competitive with CMB bispectrum measurements [4]. 
For other models of primordial non-Gaussianity however. 
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comparisons between the high-peak expectation [Eq. (3)] 
and the simulated NG bias have revealed large discrep- 
ancies in the magnitude of Abj [6, 13, 14]. 

In [8] , we calculate the scale-dependent NG bias using 
two different formulations of the peak-background split. 
The first generalizes the approach adopted in [4, 7]; the 
second utilizes conditional mass functions. Both methods 
arrive at the same expression on large scales and thus 
validate the robustness of the result. Let us briefly review 
the second approach, which is inspired by a derivation of 
the Gaussian peak bias factors [12]. The conditional mass 
function n{M, Rs\Si, Ri) gives the mean number density 
of halos of mass M and radius Rs inside a large region of 
size Ri ^ Rs and overdensity Si. The conditional mass 
function encodes information on the halo bias parameters 
[10, 11]. On taking the limit i?; — > oo and expanding the 
conditional over-abundance of halos. 



niM,R4Si,Ri) 



(5) 



in powers of the large-scale perturbation Si, we are able 
to read off the NG bias factors. In particular, the scale- 
dependent NG contribution to hf is at leading order 



Ab 



(pbs) 



(fc) 



4 J-W(fc) 
(iV-1)! MM 



(6) 



e,(fc) = bN-2Sc + bN-3 3 - AT + 



where ^at is the 7V-th order, Gaussian Lagrangian bias 
(with bo = 1). The first term {hM~2Sc) reduces to Eq. (3) 
for high mass halos since ^at ~ (v/cros)^ when v ^ 1. 
Therefore, it exactly reproduces the previous results of 
[4, 6, 7, 13]. The second term, which is derived for the 
first time in [8], only vanishes in the local, quadratic /nl 
model and on large scales. Although it becomes sub- 
dominant in the high-peak limit {bN-i/bN-a ^ 1), it is 
significant for most relevant peak heights. 

To decipher the origin of this new term, it is useful 
to recall the basic mechanism by which primordial non- 
Gaussianity affects halo clustering. In the presence of 
a primordial bispectrum, the small-scale density field is 
modulated by a long-wavclcngth (fc < 0.01/i/Mpc) per- 
turbation of the potential via [4, 7] 



(7) 



Thus, instead of being a constant throughout space, &Qg 
and hence v vary from point to point. Since the local 
abundance of halos depends on the amplitude of small- 
scale density fluctuations (it is commonly considered to 
be proportional to some multiplicity function f(y)), the 
change in a^s [Eq. (7)] contributes to the variation in the 
overdensity Sh of halos. Most importantly, this modula- 
tion of (Tqs induces a dependence of Su on (j), rather than 



Sm- This explains why the scale-dependent NG bias gen- 
erally increases towards larger scales (but note the ad- 
ditional fc-dependence brought by !Fs^\k)). As shown 
in [8], similar arguments can be made for a generic pri- 
mordial iV-point function. In this case, long-wavelength 
potential perturbations locally induce a reduced moment 
of order [N — 1) in the small-scale density field Js, which 
is given by 



5i^-i)(fc) ^ 



N-l 

'as 



= 4a 



Z-N 



J-W(fc)^Kk) . (8) 



For example, a primordial trispectrum generates a local 
skewness. Therefore, we can think of the effect of a pri- 
mordial A^-point function as locally rescaling the local 
significance according to [8] 



y 1 



r6A.-3^i^Hfc)<^,(k) 



{N - l)!"""""" " ^■-^■^'v"/ J ■ (9) 

Clearly, the modulation of the variance, skewness etc. of 
the small-scale density field not only changes the signif- 
icance V that corresponds to a mass M, but also affects 
the significance interval dv which a fixed mass bin dM is 
mapped to. Hence, long-wavelength perturbations of the 
Jacobian d\x\v/d\x\M generate a dependence of the non- 
Gaussian halo bias on the derivative of the shape fac- 
tor J-s w.r.t. Rs- This term must be present because 
halos in N-body simulations and, to a lesser extent, in 
galaxy surveys are identified by mass. In contrast, in the 
thresholding approach the two-point correlation function 
of halos above a mass threshold M is associated with that 
of Lagrangian regions above a threshold S^. in the linear 
density field smoothed on a fixed scale Rg{M). This is 
the reason why this approach cannot recover the second 
term in the expression of es{k) [Eq. (6)]. 

In the remainder of this Letter, we shall compare the 
prediction Eq. (6) to simulation results in the large-scale 
limit fc — > 0. In this regard, it is convenient to express 
both the predicted and simulated scale-dependent NG 
bias in units of the high-peak expectation Eq. (3), i.e. 
we plot Abi/Abj^^. Our prediction for this ratio is 



Ab 



(pbs) 



(fc) 



Ab 



_ , _ (10) 
.^^^\k) \<^OsJ biSc 

for a given primordial A'^-point correlation function. Note 
that, in the models we consider below, €s depends very 
weakly on k for wavenumbers k < 0.01 /i/Mpc. Hence, 
Eq. (10) predicts a mass-dependent correction to the am- 
plitude of the scale-dependent NG bias on large scales. 

We consider the following three models beyond the 
scale-independent local /nl for which simulations have 
been performed: 

• scale-dependent /nl [13] 



(ki,k2,k3) 



/NL(fcp) 



fcl 



P2P3 



(5 perm.) 

(11) 
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FIG. 1: Ratio of the large scale NG halo bias measured in 
simulations with non-Gaussianity of the scale-dependent local 
type (n/ = ±0.6) [13], to that predicted by the high-peak 
approximation. Filled and open symbols show the simulation 
results, whereas the solid and dashed curves indicate the new 
theoretical expectation [Eq. (10)] at z — and 1, respectively. 
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FIG. 2: Same as Fig. 1 but for the orthogonal bispectrum 
shape [14]. Different symbols correspond to various halo mass 
bins as indicated in the figure; colors indicate the redshifts 
z — (magenta), 0.67 (blue), 1 (green) and 1.5 (red). The 
solid curve shows the new theoretical prediction assuming a 
halo mass M = 5 x 10^^ Mq//i. 



• orthogonal bispectrum [14] 

ej,^)(ki,k2,k3) - 6/nl[PiP2 + (2 perm.) + 3(PiP2P3)'/' 
- {pI^'p^/'Ps + (5 perm.))] (12) 

• local f^NL model [6] 

4^'(ki,k2,k3,k4) = 6gf,L[PiP2P3 + (3 perm.)] , (13) 

where Pi = Ptp{ki). The shape factor for each of these 
models on large scales respectively asymptotes to 

-Ff )(fc) = /NL(fcp)fc;"^' f — V , « = nf/2 (14) 

Ff\k) = -3/nl (—\ fc-'" , a = (ns ~ 4)/6 (15) 

\O-0s / 

:F^''\k) = \9^i.^lsfl,, (16) 
where the spectral moments are defined as 

aL^E j^^k^-p,{^k)Ml{k). (17) 

In Eq. (16), S].{^^ is the skewness of the density field 
smoothed on scale Rg in a local quadratic model with 
/nl ~ 1- Using these expressions, it is straightforward 
to evaluate the ratio Eq. (10). 



Fig. 1 shows the simulations results of [13] for the lo- 
cal, scale-dependent /NL-type model [Eq. (11)] relative 
to the high-peak prediction, as a function of the Gaus- 
sian halo bias . The values of fcp = 0.04 Mpc~^ and 
fisiLikp) = 630 are fixed, while the spectral index is 
Uf = ±0.6. Clearly, the data points are inconsistent 
with the high-peak prediction in the range bf ^ 4. The 
improved theoretical prediction presented in this paper 
explains the observed deviation of the high-peak predic- 
tion for the z = case (filled symbols, solid lines) and for 
the z = 1 case and negative nj (open symbols, dashed 
lines). The only exception is for nf = 0.6 at z = 1, 
where the N-body data appears consistent with the high- 
peak prediction. This case deserves further investigation. 
Note however that the simulation measurements arc cor- 
related, since they were estimated from a common simu- 
lation volume. Hence, the significance of this departure 
is not straightforward to estimate. 

The results from simulations of the orthogonal bispec- 
trum [Eq. (12)] are shown in Fig. 2 [14]. Symbols repre- 
sent the measured NG halo bias from simulations relative 
to the high-peak prediction from three realizations with 
/nl = —250 and two with /nl = —1000. The results 
were averaged over wavenumbers with k < 0.1 /iMpc~^. 
We show our prediction (Eq. (10) with iV = 3 together 
with Eq. (15)) as the solid curve, for a median halo mass 
M = 5 X 10^"^ Mq/H. Again, the improved theory is con- 
sistent with the N-body data, which convincingly shows 
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FIG. 3: Same as Fig. 1 but for tlie local cubic model. Triangle 
and square symbols indicate results for (jnl = —10*' and 10* , 
respectively. Different colors show different redshifts spanning 
the range < z < 2. 

a strong suppression for &f < 2. While the error bars 
are significant, the strong mass- and redshift-dependence 
predicted by Eq. (6) appears to be supported by the data. 

Finally, Fig. 3 shows the ratio of the simulated large- 
scale bias of the local ^nl model [Eq. (13)] to the high- 
peak prediction. The simulations used values of ^nl = 
±10^ [6]. Error bars indicate the scatter among 5 real- 
izations. Clearly, the measured values of Abi lie far be- 
low the high-peak prediction. While they increase mono- 
tonically with 6f , Abi/Ab|''P^ never reaches unity even 
for the most biased samples. The solid curve represents 
the improved theoretical prediction, Eq. (10) with the 
shape factor Eq. (16). For this Figure, we have com- 
puted the Gaussian peak-background split biases &i and 
&2 from a Sheth-Tormen mass function [11]. The pre- 
diction matches the measurements well, though it some- 
what underestimates Abi for bf ^ 2; an improved match 
might be achieved by measuring the 62 of halos directly 
from simulations, rather than using the fitting function 
of [11]. Note also that the measurements are noticeably 
larger for ^nl — —10^, suggesting that second order con- 
tributions in are important. This is presumably due 
to the fact that the Gaussian bias factors &i and &2 re- 
ceive scale-independent NG corrections which depend on 
the nonlinear parameter ^nl- 

In summary, we have presented a new formula for the 
NG bias based on an improved peak-background split ar- 
gument [8]. We have obtained an additional term which 
can be interpreted as the effect of non-Gaussianity on the 
mapping between the significance v and the mass of halos 



M. This new term vanishes only for the local quadratic 
model with constant /nl, and is significant for all other 
models considered here. We compare our theoretical pre- 
dictions to the scale-dependent bias measured from NG 
N-body simulations, and show that it is in good over- 
all agreement with the measurements, in contrast to the 
existing formulae based on the statistics of high peaks. 
The improved theoretical predictions presented here will 
enable accurate upper limits on, or measurements of, 
general shapes of non-Gaussianity. In particular, our 
results are of relevance to the recent measurements of 
[15], who employed the high-peak prediction, and are ex- 
pected to increase their upper limits by factors of order 
unity. We expect similar mass-dependent corrections to 
the quadratic and higher-order bias factors which should 
significantly affect the halo bispectrum [8]. Accounting 
for all these corrections will be essential to obtain ac- 
curate constraints on primordial non-Gaussianity from 
large-scale structure. 

We would like to thank Sarah Shandera and Chris- 
tian Wagner for providing the data shown in Fig. 1 and 
Fig. 2, respectively. DJ and FS are supported by the 
Gordon and Betty Moore Foundation at Caltech. VD is 
supported by the Swiss National Foundation under con- 
tract No. 200021-116696/1 and FK UZH 57184001. 



[1] E. Komatsu, et al., Astrophys. J. Supp.192, 18 (2011). 

[2] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[3] V. Desjacques and U. Seljak, Classical and Quantum 

Gravity 27, 124011 (2010). 
[4] N. Dalai, O. Dore, D. Huterer, and A. Shirokov, Phys. 

Rev. D77, 123514 (2008); A. Slosar, C. Hirata, U. Seljak, 

S. Ho, and N. Padmanabhan, JCAP 8, 31 (2008). 
[5] S. Matarrese and L. Verde, Astrophys. J. Lett. 677, 

L77 (2008); L. Verde and S. Matarrese, Astrophys. J. 

Lett. 706, L91 (2009); 
[6] V. Desjacques and U. Seljak, Phys. Rev. D81, 023006 

(2010). 

[7] F. Schmidt and M. Kamionkowski, Phys. Rev. D82, 
103002 (2010). 

[8] V. Desjacques, D. Jeong, and F. Schmidt, Phys. Rev. 

D84, 063512 (2011). 
[9] S. Matarrese, L. Verde, and R. Jimenez, Astrophys. 

J. 541, 10 (2000); M. Lo Verde, A. Miller, S. Shandera, 

and L. Verde, JCAP 4, 14 (2008); M. Maggiore and A. 

Riotto, Mon. Not. R. Astron. Soc.405, 1244 (2010). 
[10] H. J. Mo and S. D. M. White, Mon. Not. R. Astron. 

Soc.282, 347 (1996). 
[11] R. K. Sheth and G. Tormen, Mon. Not. R. Astron. 

Soc.308, 119 (1999). 
[12] V. Desjacques, M. Crocce, R. Scoccimarro, and R. K. 

Sheth, Phys. Rev. D82, 103529 (2010). 
[13] S. Shandera, N. Dalai, and D. Huterer, ArXiv e-prints 

(2010). 

[14] C. Wagner and L. Verde, ArXiv e-prints (2011). 
[15] J.-Q. Xia, et al., ArXiv e-prints 1104.5015 (2011). 



